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AN  AUTOMATED  SYSTEM  FOR  CHARACTERIZING  LASER  PULSES 


I.  INTRODUCTION 

1.  Subject.  This  report  describes  a  scheme  to  calculate  laser  pulse  power,  energy,  and 
other  parameters  by  means  of  a  waveform  digitizer  interfaced  with  a  computer.  More 
specifically,  we  shall  characterize  a  system  which  analyzes  a  pulse  emanating  from  a  CO., 
Transversely  Excited  Atmospheric  (TEA)  laser.  The  pulse  is  detected  by  a  Rofin  Model  7441 
photon  drag  monitor,  and  the  waveform  is  digitized  and  processed  by  a  Biomation  Model  6500 
500  MHz-Digitizer  and  a  Hewlett  Packard  Model  9836  computer,  respectively.  Although  it  is 
this  specific  system  that  is  described,  it  is  the  authors'  hope  that  the  techniques  utilized  will 
prove  to  be  of  general  value  to  laser  users.  In  this  spirit,  we  have  included  also,  a  sample 
BASIC  program  written  for  the  HP  9836  and  shall  explain  it  in  depth. 

2.  Background.  In  order  to  produce  a  better  TEA  laser  or  to  make  performance  calcula¬ 
tions  for  existing  TEA  lasers,  it  is  necessary  to  have  an  accurate  means  of  determining  pulse  in¬ 
tensity  and  energy.  There  are  several  types  of  commercially  available  detectors  for  these  ap¬ 
plications.  We  shall  describe  how  two  of  them  (a  photon  drag  device  and  a  pyroelectric  detec¬ 
tor)  can  be  utilized  in  an  automated  data  recording  and  processing  system. 

The  photon  drag  detector  utilized  in  this  experiment  is  a  fast  ( <  1  ns  rise  time)  doped 
germanium  device  which  produces  a  voltage  proportional  to  the  instantaneous  intensity  of  a 
laser  pulse.  This  is  displayed  on  an  oscilloscope  to  obtain  a  record  of  the  temporal  history  of 
the  pulse.  A  typical  constant  of  proportionality  is  given  by  the  manufacturer*,  but  an  exact 
value  depends  on  beam  profile,  input  intensity,  and  experimental  set-up.  It  has  even  been 
observed  that  changes  in  the  internal  battery  supply  voltage  will  affect  the  constant.  Therefore, 
it  is  necessary  to  determine  the  photon  drag  constant  for  each  experimental  set-up  and 
periodically  check  this  value  to  update  it  if  necessary. 

The  pyroelectric  detector  employed  in  this  set-up  is  a  Laser  Precision  Model  314  device 
which  is  essentially  a  capacitor  formed  by  depositing  metal  electrodes  on  pyroelectric  material. 
When  the  radiation  strikes  the  coated  surface,  the  heat  produced  causes  a  polarization  change 
which  gives  rise  to  an  output  voltage  proportional  to  input  energy.  The  constant  of  propor¬ 
tionality  can  be  determined  by  a  calorimeter  (e.g.,  Scientech  Model  362),  and,  therefore,  this 
detector  can  provide  a  relatively  quick  means  of  measuring  the  laser  pulse  energy. 

Determining  the  laser  pulse  energy  and  intensity  with  these  detectors  can  be  done 
manually  for  a  single  pulse;  but  for  statistically  significant  numbers  of  pulses,  this  method 
becomes  extremely  inconvenient  at  best.  We  shall  describe  a  practical  method  of  interfacing 
the  photon  drag  and  pyroelectric  detectors  with  a  computer  to  calculate  real  time  pulse  energy 


714  kW/V  for  (hi*  Model  7441  with  a  built-in  X10  amplifier. 


and  intensity  as  well  as  other  parameters  concerning  the  temporal  history  of  the  pulse.  For  ex¬ 
ample,  a  typical  TEA  laser  pulse  has  the  intensity  profile  as  shown  in  Figure  1.  A  high  intensi¬ 
ty  spike  containing  20  to  80  percent  of  the  total  energy  if  followed  by  a  long,  low  intensity  tail. 
Representative  values  for  the  spike  half  width  (FWHM)  and  the  total  duration  of  the  pulse  are 
70  ns  and  2  /is  respectively.  In  addition  to  the  intensity  and  energy,  the  system  described 
enables  the  operator  to  measure  the  half  width,  total  duration,  and  amount  of  energy  con¬ 
tained  in  the  spike.  If  a  number  of  pulses  are  to  be  collected,  the  average  and  standard  devia¬ 
tion  of  the  ensemble  energy  and  intensity  can  be  displayed. 

II.  THEORY 

The  following  analysis  is  used  in  order  to  find  the  peak  power  and  energy  of  the  laser 
pulse.  The  peak  power  P  incident  on  the  photon  drag  detector  is  related  to  the  output  voltage 
Vpdby 


p  =  KPd  V 


(1) 


where  K^,  is  the  photon  drag  constant  (watts/volt). 

The  total  energy  of  the  pulse  E  is  related  to  the  power  by 

T 


e=7 


pdt. 


where  T  is  the  total  duration  of  the  pulse. 


Thus, 


T 

E=/  w 


(2) 


(3) 


If  we  assume  that  the  incident  intensities  are  low  enough  so  that  is  independent  of 
time,  then  the  photon  drag  constant  may  be  expressed  as 


Kpd  =  E/A, 


(4) 


where  A  is  the  area  under  the  voltage-time  curve. 


If  a  pyroelectric  detector  is  used  to  determine  the  energy  of  the  pulse,  equation  (4) 
becomes 


(5) 
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Figure  1.  Typical  TEA  fa*ar  pulta. 


where 


K  =  pyro  eon.stanl  (volts/joule). 

and  V  =  maximum  voltage  produced  by  laser  pulse  impacting  the  pyro  detector. 

The  pyro  constant  can  easily  be  determined  by  averaging  a  suitable  number  of  laser  pulses  of 
known  energy  (e.g.,  from  calorimeter  measurements)  and  bv  dividing  that  value  by  the  max¬ 
imum  voltage  produced  by  these  pulses.  The  determination  of  the  area  A  may  be  straightfor¬ 
ward  in  principle  but  in  practice  may  prove  to  be  quite  difficult  to  do  accurately.  For  example, 
triangular  approximations  for  the  pulse  are  too  crude  and  the  presence  of  a  long,  low  intensity 
tail  that  contains  a  significant  portion  of  the  energy  makes  planimeter  measurements,  also,  in¬ 
herently  inaccurate.  It  is  precisely  this  consideration  that  leads  one.  to  the  choice  of  digitizing 
the  output  of  a  photon  drag  detector  and  using  a  computer  to  numerically  integrate  under  the 
voltage-time  curve.  At  any  rate,  it  is  easy  to  see  how  the  photon  drag  constant  can  now  be  ob¬ 
tained  from  equation  (5)  and  therefore  how  the  energy  and  peak  power  can  be  measured  for 
any  subsequent  pulse  by  the  application  of  equations  (3)  and  (1)  respectively. 

II.  EXPERIMENTAL  SET-UP 

Figure  2  illustrates  the  arrangement  of  the  apparatus  schematically.  The  TEA  laser  is  of 
conventional  design  and  is  preionized  by  means  of  an  array  of  sparks.  The  Rogowski-profiled 
electrodes  are  separated  by  a  gap  of  1  cm  and  are  approximately  11  cm  long  by  3  cm  wide. 
When  this  device  is  fired  by  means  of  a  spark  gap  discharging  a  10-nF  capacitor  charged  to  20 
to  35  kV,  it  can  produce  output  energies  in  the  30  to  150  mj  range  and  peak  powers  from  100 
to  1000  kW.  Note  that  the  photon  drag  monitor  is  constructed  with  an  internal  beam  splitting 
arrangement  that  allows  part  of  the  incident  beam  to  be  transmitted  directly  to  the  pyroelectric 
detector.  This  allows  the  photon  drag  detector  to  be  conveniently  calibrated  without  moving 
the  arrangement.  As  previously  mentioned,  the  digitizer  is  a  Biomation  Model  6500  500-MHz 
device,  and  the  computer  is  a  Model  9836  manufactured  by  Hewlett  Packard.  The  digitizer  is 
simply  a  solid  state  analog-to-digital  converter  with  a  memory.  It  will  record  at  sample  rates  up 
to  2  ns/sample  and  store  1024  samples  (a  6-bit  resolution/sample).  The  interface  is  a  Biomation 
Model  4880  IEEE  488  coupler  that  converts  the  digitizer  parallel  binary  format  to  an  IEEE 
488  ASCII  format.  Implementation  of  our  relatively  long  and  slow  BASIC  program  limits 
pulse  repetition  rates  to  0.7  Hz  or  less,  but  for  most  applications,  this  is  not  a  severe  restriction. 

The  biggest  problem  with  this  set-up  is  not  directly  caused  by  the  digitizer/computer  ap¬ 
paratus  but,  instead,  by  the  device  which  it  is  measuring,  namely  the  laser.  Because  of  the 
nature  of  the  discharge  in  TEA  laser,  it  is  inherently  electrically  noisy.  For  many  commercial 
lasers  this  is  no  big  obstacle  since  ordinary  metal  enclosures  provide  sufficient  shielding. 
However,  laboratory  breadboarded  designs  are  typically  unshielded,  and  obtaining  a 
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favorable  signal-to-noise  ralio  imi\  pose  a  problem.  Shielding  a  research  laser  itself  mav  no!  he 
practical  since  circuitry  may  need  to  he  changed  often,  currents  and  voltages  must  he 
monitored,  and  electrical  discharges  often  need  to  he  visually  observed.  But  without  shielding, 
noise  levels  may  he  high  enough  to  activate  the  registers  in  the  digitizer  even  w  ith  no  cable  con¬ 
nections. 

In  order  to  deal  with  this  problem,  we  have  found  it  expedient  to  utilize  several  tech¬ 
niques.  First,  simply  moving  the  digitizer/computer  3  to  4  m  away  helj»ed  greatly.  The  extra 
cable  length  then  required  the  use  of  double-shielded  coaxial  cable.  As  a  further  precaution, 
the  box  was  constructed  of  plexiglass  and  lined  with  copper  mesh.  When  the  digitizer  and  in¬ 
terface  were  placed  inside,  the  combination  of  these  measures  provided  sufficient  shielding  for 
all  except  the  most  severe  conditions  encountered,  when  extremely  high  ( >  32  kV)  voltages 
were  needed  to  charge  capacitors  during  certain  high  pressure  experiments.  In  that  case,  the 
pulse  forming  network  was  enclosed  in  a  shielded  box.  Workable  signal-to-noise  ratios  were 
then  restored. 

IV.  COMPUTER  PROGRAM 

For  the  purpose  of  clarity  to  both  computer  programmer  and  user,  two  block  flowcharts 
are  presented.  The  first  (Figure  3)  is  a  simplified  input/output  (I/O)  flowchart  for  the  user. 
After  the  program  is  running,  the  user  should  be  able  to  follow  this  flowchart  with  verv  little 
knowledge  of  the  algorithms  used  to  determine  the  output.  The  second  flowchart  (see  Appen¬ 
dix.  Figure  A-l)  is  for  the  programmer  and  gives  a  much  more  detailed  explanation  of  the 
algorithms  used  in  the  determination  of  pulse  parameters  as  well  as  input  and  output 
algorithms.  Although  the  program  was  written  in  BASIC  for  the  HP  9836  computer  and 
Biomation  6500  digitizer,  the  calculation  algorithms  are  universal,  however.  I/O  algorithms 
will  probably  be  different  for  other  systems.  A  sample  computer  program  listing  and  output  ex¬ 
amples  are  included  in  the  Appendix. 

^  hen  the  program  is  run,  the  photon  drag  constant  and  py  roelectric  constant  are  set  to 
default  values.  These  initial  numbers  are  not  important,  because  thev  can  be  updated  after  the 
sample  of  pulses  has  been  recorded.  The  digitizer  input  voltage  and  time  base  settings  must  Ik* 
entered  into  the  program.  This  is  done  automatically  by  manually  arming  and  triggering  the 
digitizer  when  it  is  prompted  by  the  program.  The  program  is  now  ready  for  pulse  waveform 
input.  W'hen  a  sufficient  number  of  pulses  has  been  recorded,  the  user  hits  Key  0  and  the  pro¬ 
gram  goes  into  a  calculations  mode.  The  initial  data  normalization  takes  about  9  s/pulse.  After 
this,  the  individual  pulse  parameters  are  calculated  and  printed  out  on  the  CRT  sequentially. 
Finally,  the  average  and  standard  deviation  of  energy  and  maximum  intensity  for  the  entire 
sample  of  pulses  are  calculated  and  printed  out. 

The  program  is  now  idling  in  an  output  mode  and  can  be  continued  by  depressing  the 
appropriate  user  defined  key.  After  the  corresponding  key  function  is  performed,  the  program 
returns  to  the  idling  point  and  awaits  another  key  command.  The  pulses  are  referred  to  hv 
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Figur*  3.  Simplified  I/O  urer  flowchart. 


number  in  the  sequence  in  which  they  were  recorded  (i.e..  the  first  pulse  recorded  is  No.  1,  etc.). 
Thus,  w  hen  prompted  to  refer  to  a  pulse,  the  operator  should  he  entering  an  integer  less  than  or 
equal  to  the  iiuniitcr  of  pulses  in  the  sample. 

A  plot  of  peak  power  vs.  time  can  he  produced  on  the  computer  CRT  by  depressing  Key 
0.  The  axes  are  sealed  automatically  hv  the  sample  intensity.  An  identical  plot  can  he  produc¬ 
ed  on  an  external  paper  plotter  hv  depressing  Key  1  A  typical  output  is  shown  in  Figure  4. 

If  a  paper  printout  is  desired.  Key  2  is  depressed  and  all  further  type  will  he  sent  to  the 
external  printer.  If  the  operator  later  wishes  the  type  to  la*  printed  on  the  CRT  again,  depress¬ 
ing  Key  2  once  more  will  cause  all  further  output  to  he  displayed  on  the  screen. 

A  permanent  record  of  any  individual  pulse  may  lie  obtained  by  pressing  Key  3.  This 
enables  the  user  to  store  a  file  of  pulse  intensity  on  disc,  where  it  can  he  retrieved  later  for  other 
analysis. 

If  the  operator  desires  to  measure  a  new  series  of  pulses,  he  presses  Key  4.  At  this  time, 
the  digitizer  sensitivity  and  time  base  settings  may  be  changed,  if  desired.  An  important 
feature  of  this  key  is  that  the  same  values  for  the  photon  drag  and  pyroelectric  constants  will 
be  retained,  whereas  they  would  have  l»een  reset  to  the  default  values,  if  the  program  had  been 
simply  started  again. 

In  order  to  obtain  a  listing  of  the  parameters  for  any  single  pulse,  the  user  need  only 
push  Key  5.  An  example  of  a  typical  output  can  be  seen  in  Figure  5.  A  printout  of  the  pulse 
rise  time,  spike  duration,  total  duration,  spike  half  width,  spike  energy,  total  energy,  and  peak 
intensity  will  be  sent  either  to  the  CRT  or  to  an  external  printer  as  specified  by  Key  2. 

Key  6  is  used  to  retrieve  the  averaged  and  standard  deviation  of  the  energy  and  peak 
intensity  for  the  entire  collection  of  pulses.  It  is  important  to  note  that  pulses  with  peak  inten¬ 
sities  of  less  than  40  percent  of  the  peak  intensity  of  the  series  are  rejected  in  calculating  this 
average.  This  feature  prevents  the  biasing  of  statistics  by  false  triggering  or  other  random 
events.  However,  all  pulses,  “good”  or  “bad,”  are  recorded  and  the  information  about  any  one 
of  them  is  still  retrievable  even  though  some  are  excluded  in  the  averaging  process.  In  addition, 
the  program  keeps  track  of  the  number  of  “bad”  pulses  so  that,  for  example,  in  an  otherwise 
well-behaved  system,  the  numlier  of  laser  misfirings  during  a  series  of  shots  can  be  obtained. 
This  feature  has  definite  utility  for  those  performing  laser  lifetime  tests. 

The  photon  drag  constant  can  be  changed  so  that  the  average  energy  of  the  ensemble 
agrees  with  that  obtained  with  cither  the  calorimeter  or  the  pyroelectric  detector.  (This  pro¬ 
bably  will  need  to  be  done  after  the  first  run  with  the  default  constant.)  When  the  operator 
presses  Key  7,  he  will  be  asked  whether  the  energy  reading  will  be  done  with  the  calorimeter  or 
the  pyrolectric  detector.  If  the  reading  is  obtained  w  ith  calorimeter,  the  operator  is  prompted  to 
enter  it  via  the  keyboard  in  millijoules.  (f  the  pyroelectric  detector  is  to  be  used,  the  operator  is 


8 


PULSE  #  1  ANAI  V  SIS 


PULSE  RISE  TIME . 45  ns 

SPIKE  HALF-WIDTH . 55  ns 

TAIL  DURATION .  245  ns  AFTER  START  OF  PULSE 

TOTAL  DURATION .  3205  ns 

SPIKE  ENERGY . 15.69  mj 

TOTAL  ENERGY .  54.99  mj 

PEAK  INTENSITY . 160.7  kW 

PULSE  #2  ANALYSIS 

PULSE  RISE  TIME . 45  ns 

SPIKE  HALF-WIDTH . 50  ns 

TAIL  BEGINS . 170  ns  AFTER  START  OF  PULSE 

TOTAL  DURATION .  3235  ns 

SPIKE  ENERGY . 12.84  mj 

TOTAL  ENERGY .  55.24  mj 

PEAK  INTENSITY . 175.3  kW 

PULSE  #3  ANALYSIS 

PULSE  RISE  TIME . 50  ns 

SPIKE  HALF- WIDTH . 50  ns 

TAIL  BEGINS .  265  ns  AFTER  START  OF  PULSE 

TOTAL  DURATION .  3205  ns 

SPIKE  ENERGY . 16.1  mj 

TOTAL  ENERGY . 54.16  mj 

PEAK  INTENSITY . 157.1  kW 

PULSE  #4  ANALYSIS 

PULSE  RISE  TIME . 50  ns 

SPIKE  HALF-WIDTH . 55  ns 

TAIL  BEGINS . 115  ns  AFTER  START  OF  PULSE 

TOTAL  DURATION . 3160  ns 

SPIKE  ENERGY . 10.75  mj 

TOTAL  ENERGY .  55,07  mj 

PEAK  INTENSITY . 175.3  kW 


STATISTICS  EXCLUDING  BAD  PULSES 

AVERAGE  MAXIMUM  INTENSITY  = . 166.9  kW 

STANDARD  DEVIATION  OF  INTENSITY  = . 7.9  kW 

AVERAGE  ENERGY  = . 55.0  mj 

STANDARD  DEVIATION  OF  TOTAL  ENERGY  = . 5  mj 

PHOTON  CONSTANT =116.88  PYRO  CONSTANT =459.13* 

TOTAL  OF  0  BAD  PULSES  OUT  OF  10  RECORDED 


Figura  6.  PuIm  parameter  listing. 


10 


directed  to  select  that  detector  via  a  SPDT  switch  (see  Figure  2)  and  to  fire  a  number  of  pulses 
into  it.  The  average  value  of  the  peak  voltage  is  automatically  calculated  and  used  to  deter¬ 
mine  the  energy  of  the  laser,  whereupon  this  information  is  used  to  calculate  the  new  photon 
drag  constant.  In  either  case,  all  future  pulse  data  will  automatically  us*  this  updated  con¬ 
stant. 


In  order  to  use  the  pyroelectric  detector,  it  is  necessary  to  calibrate  it.  Pressing  Key  8  will 
facilitate  this  task.  The  user  is  directed  to  let  several  pulses  impinge  on  the  detector,  and  as 
previouslv  stated,  the  output  is  averaged.  Then,  a  calorimeter  is  used  to  find  the  actual  energy 
and  this  value  is  entered  via  a  prompt.  The  new  pyro  constant  is  displayed  and  used  thereafter. 

As  indicated  in  the  flowchart  (Figure  3),  Key  9  is  used  to  end  the  program.  The  user 
should  be  aware  that  this  key  must  not  be  used  to  simply  enter  a  new  series  of  shots  (Key  4 
should  be  chosen  for  that).  Rather,  this  key  reloads  the  internal  disc  program  directory  and, 
thereby,  facilitates  the  easy  transfer  to  another  program. 

V.  RESULTS 

As  noted  earlier,  we  have  used  this  instrumentation  to  measure  pulse  characteristcs  over 
a  wide  range  of  intensities  and  energies.  Generally  speaking,  the  results  have  been  reasonably 
self-consistent  and  useful.  As  an  indication  of  the  former,  consider  the  following  results.  The 
output  from  a  Marconi  T250  TEA  laser  was  measured  with  the  photon  drag/digitizer  combina¬ 
tion.  Usually,  great  pains  are  taken  to  ensure  optical  alignment  and  rigidity  of  all  components. 
For  example,  an  optical  table  and  rail  are  used  together  with  a  clamping  arrangement  that 
prevents  relative  motion  between  the  laser  under  test  and  the  detector.  This  time,  however,  we 
deliberately  chose  a  potentially  more  unstable  configuration  to  test  the  sensitivity  of  this  re¬ 
quirement.  The  laser  was  set  on  a  lab  jack,  and  the  detector  was  mounted  on  a  simple  holder. 
Neither  device  was  fixed  to  the  table.  Four  separate  measurements  of  ten  shots  each  were  made. 
After  each  measurement  the  detector  was  removed  and  replaced  before  the  next  set.  The  posi¬ 
tioning  was  aided  hv  means  of  a  carbon  block  that  sparked  when  struck  by  a  pulse.  No  attempt 
was  made  to  keep  the  distance  from  the  laser  to  the  detector  constant.  The  measured  intensities 
are  given  in  Table  1. 

Table  1.  Peak  Powers  Measured  in  Four  Separate  Runs  of  Ten  Shots  Each 

Trial  No.  Power  (kW) 

1  176.2  ±  5.4 

2  172.9  ±  8.4 

3  168.7  ±13.2 

4  166.1  ±11.4 


H 


Since  the  averages  are  themselves  only  spread  approximately  3  percent  about  a  central 
of  171  kW,  the  reproducibility  of  this  method  of  measurement  has  been  demonstrated. 

It  should  be  noted  that  the  first  three  trials  were  taken  with  a  digitiser  sampling  rate  of  3 
ns  channel.  This  relatively  crude  measurement  was  necessary  in  order  to  make  sure  that  the  en¬ 
tire  pulse  would  be  recorded.  That  is,  since  the  pulse  lasted  almost  3400  ns,  then  the  fastest  set¬ 
ting  of  2  ns/channel  X  1024  channels  would  not  be  sufficient  and  the  next  setting  (5  ns)  had  to 
be  used.  If  one  chooses,  however,  to  use  the  faster  setting,  no  useful  information  is  lost  except 
for  the  total  energy  and,  of  course,  the  total  duration.  (Trial  No.  4  uses  the  2  ns  setting.)  On  the 
other  hand,  even  using  the  more  coarse  setting  of  5  ns  does  not  drastically  change  the  time 
measurements.  For  example,  using  the  2-ns  setting  gave  a  halfwidth  measurement  of  50  ±  6 
ns  while  the  5-ns  setting  yielded  53  ±  6  ns  for  the  same  measurements.  The  results  using  the 
digitizer  compared  favorably  with  those  obtained  by  analyzing  oscilloscope  traces;  e.g.,  see 
Table  2. 


Table  2.  The  Halfwidths  and  Rise  Times  of  Pulses  as  Measured 
With  a  Digitizer  and  an  Oscilloscope 


Time 

Digitizer 

Oscilloscope 

T^  (ns) 

48  ±6 

45  ±  4 

Tfise  (ns) 

47  ±3 

43  ±1 

The  oscilloscope  used  had  a  quoted  bandwidth  of  400  MHz  and  was  a  Tektronix  7834 
storage  device  with  a  7A19  amplifier  and  a  7B85  time  base.  The  measurements  were  made  us¬ 
ing  the  10-ns/div  sweep.  Digitizer  measurements  were  taken  of  40  pulses,  while  six  traces  were 
used  for  the  oscilloscope  data.  These  results  indicate  that  the  digitizer  can  consistently  produce 
measurements  that  agree  with  those  obtained  in  a  more  conventional  manner  to  within  10  per¬ 
cent. 
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VI.  CONCLUSIONS 


It  has  been  demonstrated  that  the  combination  of  a  digitizer  and  computer  is  useful  in 
analyzing  the  Qutput  of  laser  pulses.  More  specifically,  we  have  shown  that  important  pulse 
parameters  including  energy,  power,  half  width,  rise  time  and  total  duration  can  be  measured 
with  reasonable  accuracy  and  consistency.  The  purpose  in  doing  this  is  not  to  suggest  that  this 
instrumentation  be  used  in  place  of  other  simpler  and  more  conventional  measuring  devices, 
but  instead  to  show  that  for  some  applications  this  technology  offers  a  fairly  straightforward 
means  of  automatically  obtaining  and  processing  large  quantities  of  data. 

Although  this  report  specifically  deals  with  the  Biomation  6500/Hewlett  Packard  9836 
combination  of  digitizer  and  computer,  it  is  hoped  that  it  will  be  at  least  of  some  value  to  those 
using  different  apparatus.  Even  in  that  case,  a  more  detailed  examination  of  the  computer  pro* 
gram  might  yield  useful  information  to  the  reader.  For  that  purpose,  we  have  included  our 
BASIC  program  listing  and  a  short  discussion  appears  in  the  Appendix. 
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APPENDIX 


This  Appendix  is  a  more  detailed  account  of  the  computer  program  including  a  com¬ 
plete  flowchart  (see  Figure  A-l)  and  a  program  listing  (Figure  A-2).  It  is  hoped  that  this  will  he 
useful  to  not  only  those  with  the  same  digitizer/computer  combination,  but  to  others  who  might 
benefit  from  exposure  to  the  algorithms  for  event  recognition,  energy  and  intensity  calcula¬ 
tions,  etc.  The  blocks  in  the  flowchart  are  numbered  sequentially,  and  further  descriptions  are 
keyed  to  these  numbers. 

Descriptions  of  the  flow  chart  blocks  follow.  The  numbers  refer  to  the  blocks,  line  numbers 
refer  to  the  program  listing,  and  variable  names  appear  in  parentheses. 

1.  Set  Constants  (lines  70  to  80).  The  program  sets  the  values  of  the  photon  drag  con¬ 
stant  (Constant)  in  kW/V  and  the  pyroelectric  constant  (Pyro-const)  in  mJ/V.  These  values  are 
typical  but  have  to  be  changed  for  each  set-up.  They  can  be  changed  after  pulse  waveforms  are 
input,  so  default  values  are  not  important. 

2.  Read  and  Decode  Digitizer  Control  Settings  (lines  140  to  610). 

a.  Format.  Data  input  from  the  6500  digitizer  is  in  the  form  of  sample  elements. 
There  are  1024  time  sample  elements  in  a  full  screen  trace,  each  with  a  corresponding  voltage 
to  create  a  stepped  waveform.  The  value  of  each  of  the  1024  samples  can  be  set  from  2  ns  to  1  s 
giving  a  full  screen  range  of  2048  ns  to  1024  s. 

There  are  252  vertical  steps  (253  positions)  possible  for  each  time  sample  ( —  128  to 
124).  The  value  of  zero  is  always  center  screen,  the  top  is  always  1 24.  and  the  bottom  is  always 
— 128.  Each  vertical  increment  is  equivalent  to  4  steps.  The  voltage  for  each  vertical  incre¬ 
ment  is  determined  by  the  input  range  setting  (Vrange).  Since  +  Vrange  corresponds  to  a  value 
of  124  steps,  each  increment  is  (4  steps/124  steps)*Vrange  =  .03  Vrange.  The  value  of  Vrange 
can  be  set  from  .25V  to  5V. 

There  is  an  offset  which  moves  the  O  V  line  within  the  field  of  —  128  to  124.  This 
offset  is  set  in  terms  of  a  fraction  of  the  input  range  (from  —  .99  to  +  .99  x  Vrange).  The  offset 
progresses  in  increments  of  4  steps  (i.e.,  —  .99  to  .99  x  124  approximated  to  the  nearest  multi¬ 
ple  of  4).  This  is  the  point  of  the  0  V  line.  An  example  (Figure  A-3)  is  shown  for  clarity. 

b.  Algorithms.  The  digitizer  control  settings  are  encoded  in  the  last  three  bvtes  of 
data  read  from  the  digitizer  (1024  to  1026).  The  code  definitions  are  explained  in  the  “Gould 
Biomation  4880  Interface  Operating  and  Service  Manual."  page  32.  and  are  reproduced  for 
reference  (Figure  A-4). 
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Figure  A-1.  Complete  block  flowchart. 
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A-2.  Program  listing. 

_  -  r.  ^  *  _■  i_  1.  '  J  h>  •_  f  N  J.  ';j  r  r\  f  .1  (l  h*  f  }  f*t 

» 

-  1  T'l  I  J  PRObWAN  v_Ai_C,Ji_A  ■;  r.'1  iff  ,  ±  i  '*  x'  Ij  i  *R.L  H  n  n,  h '  .  tLn  _Ui*  vb. 

!  ANI\GTWER  :".-LSE  -ARAMETERS  FROM  A  V  .,<=  i  t  wAVEFGRM  FROM  A  PHOTON 

jU  ?  DRAG  Dc.it.bfUR  AND  HN  tUitRoY  REHt  a.M*o  rRiiN  A  LfiL'iRlMETc.n  .  PR'  ■  ‘  -  a  r-  *■ 

!  INTEGRATES  UNDER  INTENSITY  vs  T I  HE  TP,  DETERMINE  ENERGY 

50  ! _ _ _ 

t'Q  * ***■•*•*  |)  CONSTANT  DtFIMITIONS  AND  SETiINGS  FROM  DIGITIZER 

70  Cons tant  =  999. AE  {CONVERSION  FACTOR  FOR  PHOTON  DRAG  < KW/VOLT ) 

80  Pyro_const=459. ’3  ’CONVERSION  FACTOR  FOR  PYROELECTRIC  <r,J/VOL7> 

51  DISP 

90  Start:  GOSUB  Offjceys  !  CLEARS  KEY  FUNCTIONS 

100  PRINTER  IS  1  !  SETS  PRINTER  INTERNAL 

MO  P r inter*  ! 

120  PRINT  CHR$(12>  !  CLEARS  SCREEN 

130  GCLEAR  !  CLEARS  GRAPHICS  1 

’  AO  GOSUB  Settings 

150  GOTO  A60 

I  GO  Sett ings: PRINT  "MANUALLY  ARM  AND  TRIGGER  THE  6500  TO  RECORD  PANEL  SETTINGS." 

170  PRINT  "(NO  SIGNAL  AT  INPUT)  THEN  <CONTINUE>" 

180  BEEP  80.. 1 

190  PAUSE 

200  PRINT  CHR$(12>  ’CLEARS  SCREEN 

210  OUTPUT  70A :  1 1  ’PUTS  6500  IN  OUTPUT  MODE 

220  ENTER  70A  USING  ’LOADS  BUFFER  ARRAY  (UNFORMATTED) 

230  OUTPUT  70A;0  {CLEARS  6500  OF  OUTPUT  MODE 

2A0  SEND  7 : UNL  ’UNLISTEN  COMMAND  TO  ALL  UNITS 

250  Irnt=B<1024>  ’INPUT  RANGE  NEGATIVE  TRUE 

260  Opt=6< 1 025)  ’OFFSET  POSITIVE  TRUE 

270  Sint=B< 1 026)  {SAMPLE  INTERVAL  NEGATIVE  TRUE 

280  Weighmg=BINEOR<BINAND<Sint.  I5>.  15)  {MAKES  PQS  TRUE  WORD  OF  NEG ‘TRUE 

BYTE 

290  RESTORE  300  ’RESTORES  DATA  (LINE  2A0> 

300  DATA  2.5.10.20.50,100.200.500.1000  !  DATA  FOR  WEIGHING  KEY 

3 1  Ov  FOR  I»0  TO  Weighing 

320  ON  ERROR  GOTO  210 

330  READ  Interval  ’DETERMINES  WEIGHING  OH  SAMPLE  INTERVAL 

340  OFF  ERROR 

350  NEXT  I 

360  Off  set  =  8lNAND(Qpt ,  1 27) «< -  1  >  '  ( BI T(Oot . ?  >  > / 1 00  {DETERMINES  FRACTIONAL 

OFFSET  FROM  BINARY 

370  FOR  1=2  TO  5 

380  IF  NOT  (BIT ( Irnt , I ) )  THEN  400  ’  BRANCHES  OUT  WHEN  BIT  IS  LOW 

3 90  NEXT  I 

*00  Range*. 25*2  (6-1)  !  VOLTAGE  RANGE  GENERATION 

410  IF  1=2  THEN  Range=5  {EXCEPTION  TO  RANGE  GENERATION 

•420  PRINT  "RANGE*"  : Range : "VOLTS" 

430  PRINT  "OFFSET*" : Of  f set  :"X  INPUT  RANGE ( "  -.Off se t -Range  : "  VOLTS)" 

*40  PRINT  "INTERVAL*”; Interval :"nS" 

450  RETURN  ! 


4b0  P  r  i  nter  =  1  ’PROGRAM  FLAG  TO  KEEP  TRACK  OF  PRINTING  DEVICE 

--70  PRINTER  IS  1  ’SETS  PRINTER  INTERNAL 

430  PRINT  "THIS  PROGRAM  WILL  ARM  THE  6500  DIGITIZER  AND  READ  THE  PULSE  DA'A" 
*80  PRINT  "AFTER  A  TRIGGER  IS  RECEIVED.  MAIN  MEMORY  HANDLES  UP  TO  120  PULSE'..’ 
,  Ch'RS  ( '  D  ■!  > 

5  00  PRINT  "PYROELECTRIC  CONSTANT*" :Py ro_const : "nj/VOLT" 

5)0  PRINT  "PHOTON  DRAG  CONSTANT*" ; Constant ; "KN/VOLT  ";CKRS(I23> 


j  Reproductd  from 
1  ban  avilabla  copy. 


Figure  A-2.  Program  listing  (continued). 


'  REALTIME  SECONDS  41 


-r\  1  r  -  ujiior .  ' 


-0  ^  =  ABS 1 C  ’  • s*  t ■*  •  28 )  {CALCULATION  CF  aP^ROX I MA 1 E  L'Fj-  SET  ■  j  hr  K. 

50  l'c_ot  ;ii‘  t  =  INT  <  R/4  )  *£  ?  IF  THE  DIGITIZER  HAS  SEEN  TRIGGER::.! 

r  <j  Of :  iec\0  THt.N  Dc_ot  f ; e 1 * -Dc_of *  -i e t  !  CHANGES  uAlt  :  ^U’l  Arts 

"  "  i  1 

'SO  IF  ri<50i)>>-12S  THEN  3<50Q  )  »B<  500  ) -256  {CONVERTS  BYTE  T0  WORD  Ir  N£G. 

c'S'ir;T  "DC  OFFSET-"  ;Dc_offset  ."PIXEL  OFFSE  ■  =  ":  l(  500  )  :  “(  S'-OU  -D  EE  WlTfi’S 

00  i~  h ELS <  Or _o f  f  se  t ~ B <  50 0 )  > > h  THEN  MO  '  IF  OFF  SETS  DUN  ’  T  .-i.-REt  r:E~  i  RIGGER 

:0  Dr _i"i f  f  s«r  r  =  9 <  500 >  !  pixel  OFFSET  is  COPREC  !  DC  OFFSET  (USUALLY  > 

20  INPUT  “ENTER  OUTPUT  OPTION  MALL  PARAMETERS  2) ENERGY  AND  INTENSITY  ONLY". 


DISP  "TURN  ON  DETECTOR  AND  <CONTINUE>  TO  TAKE  DATA" 
PAUSE 


t  -mf  **■!**•  ■!*  !+•#■*■■**•*  * 


ijf  PHRAMhTER  I NPu '*■*•*■  *.»*•• 
_  BEGINING  OF  PULSE  INPUT  r-ijDE 


ON  '< E 7  0  LABEL  "END  INPUT"  GOTO  Calc  !  ENDS  PULSE  INPUT  MODE 
PRINT  CARS  < ; 2  >  !  CLEARS  SCREEN 

PRINT  »«iCLl  IN  INPUT  MODE...  READY  FOR  PULSE  INPUT" 
lain :  1  THIS  SECTION  USED  FQP  MAIL.  MEMORY  STijRAbF  ARRAY 

PRINT  .M  ING  .  -K"  : CHRS<  132),  "HIT  <0  TO  END  PULSE  INPU’  MODE  .  "  .  S'HRS<  1 28) 
INTEGER  A <  120. 1  02*3 >  ,B<  1 026 >  .Parameter  <6. 1  20  >  .  BaU_puise <  ’  20  ) 

DIM  120)  . AspiLe(  120)  .Tenpyt  1024) 

Pulse" l  !  NUMBER  OF  PULSES  RECORDED 

GUTPLT  '0-:o  1  ARMS  THE  6500  DIGITIZER 

OUTPUT  !  CLEARS  ARM  COMMAND 

OUTPUT  *  PUT  6500  INTO  OUTPUT  MODE  (SINGLE  DIGIT  FORMAT ) 

DISABLE  !  DISABLE  USER  DEFINED  KEVS 

ENTER  704  USING  !  ENTER  ARRAY  B<1024)  UNFORMATTED  (FAST) 

OUTPUT  704:0  !  CLEARS  THE  6500  OUT  OF  DIGITAL  OUTPUT 

Pu.Ise”Pu!se+ 1  !  TALLY  OF  PULSES  ENTERED 

DISP  Pulse:"  PULSES  STORED  " 

FOR  1=1  TO  1023 

A (Pul ;e . I >  =B< I >  !  MOVE  BUFFER  ARRAY  INTO  MAIN  ARRAY 

NEXT  I 

ENABLE  !  ENABLE  USER  DEFINED  KEYS 

IF  Ru. !■=.-><  1  20  THEM  750  !  MAIM  MEMORY  CAM  QMl Y  HANDLE  120  PU-.2ES 

PRINT  "MO  '-'ORE  ROOM  IN  MEMORY.  PRESS  CONTINUE  FOR  CALCULATIONS" 

PAUSE!  END  OF  PULSE  INPUT  MODE  _ _  _  _ 


:  Calc  :!!!!'  r  1  '!!!!!!!!!  !  CALCULATIONS  FOR  PULSE  PARAMETER.;.  !  !  !  !  1  1  1  !'!:!!!!!!  ' 

! . OFF  ’  KEY  0 . 

1  PRINT  CHRS(II')  ! CLEAR  SCREEN 

i  PRINT  "IN  CALCULATION  MODE  . " 

1  PRINT  "NOW  PROCESSING  ": Pulse:"  PULSES.” 

PRINT  USING  "-K":CHRS< 130). "PLEASE  WAIT  » ! ! "  .  CHR5 < 1 20 > 

DISP  "  " 

Tcal c  =  T IMEDA  i  E  ! REAL  TIME  IN  >tCl‘NDS 

SEND  7  :  IJNL  {COMMAND  UMLISTEN  TO  ALL  UNITS.  RESE  T  GFI2  A :  :  N  :  M  .M 
Maxval =0  ‘MAXIMUM  VALUE  OF  ALL  PULSES  RECORDED 

J  "  Jr  -  ,  Pulse  !  LUO*’  FUR  El  ft-  - "  F'IJl  :■  : 


C  1023 

(=127  THEN  1040 
'  > -256 

'  - < I , J) -Dc_cf f  set ) 
‘a.tvai  HEN  1070 


!  ,jl!  p  ~  OR  EACH  PIXEL  OF  E  •'.<>! 

? CONVERT  TO  WORD  2  '  S  (0^..- 'ML 
'  f'R '”M  BYTE  2 ' 2  .  .v;  .  ■, 

1  A:  Mi  .  -OR  DC  0  -  ••  ' 

•  -ND  RECTTf  Ir  s  -''UL'.'E 
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u  i.! 

;■  s !  i 


j  1 1 
n  ' ; 

:  1 70 
’  I  80 
'  190 
1  200 


.0 


290 

2b  0 
270 
2 3  0 
230 
300 
2 1  0 
320 
330 

?  40 


.■ax  .  a :  -  1 
■|£.<*  j 
NEXT  I 

»  ,  ^  #  ^  +  +  ^ 

Ba-.-O 
A  v  n  a  x  -  0 
F  OP  =■»’  ' 
A  ci.X  ~  ‘ 1 
T s i ar  t  *  0 


Figure  A-2.  Program  listing  (continued). 

■  I(:L*S  T-:£  *■!»->: I'*- .  .  'EM 


;Lv';.L'i.’  ji-  :.c.r"~'n  -  ->  ihl  ‘  -  :■  - ■* 

HUMBER  M  csAl*  P J L cJ;<_..iUN  •  lK"  j 
AVERAGE  MAXIMUM  INTENSITY 

!MAIN  LOOP  3 7 ART  FOR  EACr.  ?.„_3E 
MAXIMUM  OF  CURRENT  PULSE 
TIME  START  qf  PUlS-  OCCURS 


t  jn  i  =  :j 

IF  (I-Ts 
IP  A<P.I 

Tatar  t - 1 
IF  A i P  ,  I 
Max  =A ( P , 
T  riax  =  I 


•  500  !  SAVES  TIME  HOT  TO  lGOF  TO  102- 

!  COULD  TRIGGER  ON  RISING  SLOPE.  BUT  NOISE  MAY  INHIBIT  TRIGGER 


tart) >200  THEN  1210 
)>0  THEN  ’210 

XMax  THEN  1240 
I> 


HEX ' 


!  IF  DATA  NOT  0  FOR  200  'HEli  PULSE  STARTED 
’LOCATE  S' ART  OF  PULSE  (TRIGGERS  OFF  0) 

’TIME  START  OCCURED  (SAMPl.ES) 

!FIMD  MAXIMUM  INTENSITY 

{MAXIMUM  INTENSITY  OF  .CURRENT  PULSE  (PIXELS) 
!  T'ME  MAX  2-CCcRED 

IF  Max > . 4*Max vai  THEN  1300  ! REJECTS  PULSE  IF  UNDER  ML;  OF  SAMPLE  MAX 

Bad = Bad V  ’COUNTS  NUMBER  OF  BAD  PULSES  ENCON7ERED 

Bad_c«uise(P  >  =  *  !  ARRAY  MARKS  BAD  PULSES  SO  NON '  T  BE  IJCL.DED  IN  STATISTICS 
PRINT  "PULSE  -".P,"  IS  A  BAD  PULSE" 

GCsT0  I  SOU  ’PROCESS  NEXT  PULSE 

IF  Outoo t  =2  THEN  1  AGO  ’  OUTPUT  OPTION  2  (ENERGY  -M  INTENSITY  ONLY) 

ledge =0  ’  INITIALIZE  TAIL  START  'I"! 

FOR  I-Tmax+20  TO  800  !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<-(  <<<<<<<<<<<<<<<<<< 
iF  A ( P  .  I -  1  !  >  > A ( P  .  I +  I  0 )  THEN  1350  !  I-  -'uLSE  RISES  - :  .  AIL  STARTED 

IE  <A(P,  T-20>~A(P, I » )>=8  THEN  1  370  ! LOCATE  SLOPE  ChA’ Gi  "OR  TAIL 

TIME  OF  TAIL  S’AP' 

TRIGGERS  ON  FlftSr  CCC'JRANCE  OF  SLOPE 
*<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
ax  OR  Tmax>Tedge  OR  Tstar t>T*age  THEN  l2bO!REjECT$  ON  TIMING 
1  INITIALIZE  END  TIME  OF  D|JLSE 

TO  1023  .'////////////////////////// ////7// //////// ////////// 

THEN  1440  ! FINDS  END  OF  PULSE  (DROP  TO  ZERO) 

•TIME  END  OCCURS 

!  TRIGGERS  ON  FIRST  ZERO  GCCURANCE 
!///// //////////////// //////////// -7'// ////////////// 
e  THEN  1260  '  "EJECTS  ON  TIMING  ERROR 

!  INI  i  IAL J.ZE  ENERGY  Ar\.-.-'. 

CALCULATE  S'-'C-F  -;,D  =>ULSE  AREAS  AND  HALF  HI D'H 


SPEEDS  UP  PROGRAM  (ENERGY  AND  INTENSITY  ONLY) 
AT  ION  AT  BEGINNING  OF  TAIL 
HE  SPIFF 

!  CONDITIONS  FOR  R 


'  350 

T  edge=  I  -  1 0 

1  MO 

GOTO  1380 

■  270 

NEXT  I 

:  330 

IF  Tstar t> 

’  390 

Tend- 1023 

■  4  o  n 

FOR  I=Teog 

'410 

IF  A<  P . I ) < 

1420 

i  erid=  I 

'430 

GOTO  ’450 

'  440 

NEXT  I 

-•50 

I  r  i  e nd  <  i  e1 

■  ago 

firea< P ) =0 

q  7  0 

i- UR  i*  istar t  1 U  ! eod- 1  ! ;  ..  :•  :• ;•  :•  :  > 

4  B  0 

A  r ea  <  P ) -  A rea  <  ? > +  A ( p , I ) *  Interval 

? USES  RECTANGLE  METHOD  Tij  DEM 

4  9  0 

IF  Outopt'2  THEN 

1560  !  SPEEDS  UP 

500 

IE  I>Tedqe  THEN  1 

520  !  ENDS  IN TEC 

510 

Mv-P  i  i  p  ( P )  -  A  r  t*  c?  ( P  ) 

!  ENERGY  IN 

*7,  ‘  i  S 

IF  I>Tnax  OR  A<P. 

I>>.5*Max  THEN  ’S-i 

S  v 

T - 1 

! RISING  Hr 

S.'io 

C  I  <  '  riax  UR  A<P. 

I  X.  5-Max  THEN  15V 

-.50 

T2»  I 

!  TRA  II  i ?»'.* 

“ThTl 

NEXT  I 

»  \  •  \  •.  •  •.  \  ■.  \  \  •  • 

:ond] 


3 ICING  HALF  MAX 
’RAILING  HALF  MAX 


'  '•  Zl.i  GO  SUB  Pr  ;  ntout 

'S80  CO  TO  '730 

-SCO  P r  i  .  n t u r :  •  »••<•  +  •(-  +  +  +  •*■  ••  •*•*  +  +  •■•  +  +  • 
1  FOn  PRINT  CHRS( 1 32) ."PULSE  ;  • 
’CIO  M  Jut  opt  =2  THE  N  )  6.7:1 

7.2 il  PRINT  "VJLSE  RISE  TIME _ 

i’R'if,  ■  "SP I».f  HAL i-  -N  In  ’  H  .  ,  . 


1  _  i  .. 

ri !  in. 


f  t  r  ►  r  • 


M-  +  +  +  4+  +  +  -‘-+*- 


Copy  available  to  D7IC  does  not 
permit  fully  legible  xeproducrUon 


Figure  A  2.  Program  lifting  (continued). 


--IN’ 

*  *  O  Z-'  ~ 

_■  L  Rrn  .  i.  . 

-  f  :  n  * 

"  ■  -  .  n  . 

ij  i’’  M  l 

IN* 

-j  t*’  -  ■'  r_ 

c'lt  SGT 

I‘r 

“  "7  0 7  -L 

energy 

-  f  - 

J  ROUND 

( Max*'.o; 

“  rt.) 

PRINT 

"PEAK 

If:TEt,:,I 

: :  o 

s  E  71 J  R ! « 

f . 

!+  +  +  <■ 

r+tt+fi' 

■730 

Bad_pu 

i  s  e  <  P  ) 

=  0 

*  7 40 

Parame 

t  e  r  (  1  , 

P )  =  T  s  t  a 

•  750 

Parameter  < 2 . 

P ) =  T max 

:  760 

Parame 

t  e  r  ( 3  . 

P  >  =  Teogt 

■  770 

Parame 

ter<  4  , 

P ) = i end 

1  780 

Parame 

ter  (5 . 

P  >  =  T 1 

'  790 

Parameter  <  6 . P)  =  T2 

1  800 

NEXT  P ! ***** 

*■*■*■*•■■■*-** 

*  st  •*  -**■*■* - 

>  ’*:qe-  ,;  tart  >*:ritef'/di  :  no 
(  *  ~ n c  'start!  “Interna-  :  ‘ ‘ n $ " 

L t\ G ij t « C<  <  (He  ik  e  <  P )  ’Cons tan t “Range/  !  24000 


vge/ 


:  DRCoND <Are«i<r>  “Const  ant  “Range/’ 24000 . 4 

3 .  4 ) 


,  4  >  : 
) :  "n 


owe : 


i.H  ‘ 


►■•••■r  +  +  -t-r  +  +f  +  +  ^^-f-rt-TT-f4Tfx  +  +  +  +  -f+-t-+  +  +  -fT-»-*f+  +  +  T  +  -*--»' 

STORES  SIGNIFICANT  TIMES  -OR  EACH  PULSE  : : : : 

•  GOOD  PULSE  IF  THIS  FAR 

t  !  TIME  START  OF  PULSE  OCCURED 

!  TIME  MAXIMUM  INTENSITY  OCCURED 
?  TIME  TAIL  BEGINS  (SPIKE  ENDS > 

'  TIME  P'JLSE  ENDS 
'  TIME  OF  RISING  HALF  INTENSITY 

•  TIME  OF  TRAILING  HALF  INTENSITY 

OF  SEPARATE  PUlSE  ANALYSIS 


1810 
:  820 
1830 


Avarea-0.  !  INITIALIZE  AVERAGE  AREA 

Avmax=0 .  !  INITIALIZE  AVERAGE  MAXIMUM 

FOR  1=1  TO  Pulse I®®®®®®®®  DETERMINE  STATISTICS  OF  ALL  C'UL$ES 


v  :ji  ci  Lai  (it  at  ,o)  aj  at  ®i  iaj  ®  ;5)  (5)  Oj  ®  ioj  9  .oi  tcu  t®  0)  to,  tot  tc;  ®  toj  m  at  ®  let  L  ®  ®  ;5j  taj  !oj  ij  toj  15)  to]  t®  ;aj  ® :®  a  ®  1  ®  §1 .5  (9)  (5  §1 1®  1®  (®  ®  (S  ®  (Si'S®  ®  ®  1®  @  1*  ®  ®  w  ®  4  ®  fi>  ®  i®  M  !®  ® 


1340 
1 850 
1 860 
•870 
•880 
1890 
’  300 
1910 
1920 
1  930 
1940 
1 950 
1  960 
'970 


•  AVERAGE  MAX  INTENSITY 
?  AVERAGE  TOTAL  ENERGY 


IF  Bad_pulse( I ) = 1  THEN  1870 

Avmax=Avmax+A(  I ,  Parameter  ( 2 , 1  > )  /  <  Pal se~Bad  > 

Avarea*Avarea+Area( I ) /<  Pulse -Bad  > 

NEXT  I 

S2=0  UNITIALIZE  SD‘2  OF  INTENSITY 

S2area=0  ! INITIALIZE  SD'2  OF  ENERGY 

FOR  1=1  TO  Pulse 

IF  Bad_pul se( I >  =  1  THEN  1940  !  DON'T  INCLUDE  BAD  PULSES  I A  STATIST'S 

S2= < A ( I , Parame ter < 2 , I ) ) -Avmax )  2/(Pulse-Bad) +S2  !  SD  2  OF  INTENSITY 

S2area= ( Area( I ) -Avarea)  2/(Pulse-Bad)+S2area  !  SD ‘2  OF  ENERGY 
NEXT  I 

S=DR0UND(S0R(S2) .4)  !  STANDARD  DEVIATION  OF  INTENSTY  (ROUNDED) 

Sarea-DR0UND(S0R(S2area>.5)  !  STANDARD  DEVIATION  OF  ENERGY  (ROUNDED) . 

*  ®  4  4  4  ®  ®  4 4 f ji 4 4 4 ® ® ® 4 •'a 4 4 to)  END  OF  STATISTICS  CALCULATIONS 


•980 

Statistics : 

PRINT  CAR'S ( 132). 

"STATISTICS  EXC 

LUDING  BAD  PULSES" .CHRS ( 1 28 > 

*  9  9  0 

G0SUB  Of '  1.  eys 

_  o  0  0 

-’RINT 
8,5) : "KU1 
PRINT 

"AVERAGE  MAXIMUM  INTEN 

SI  T  Y= . 

. . " : DROUND ( Avnax*Cons tan  t * Rang 

j  0  1  0 

"STANDARD  DEVIATION  OF 

INTENSITY- _ 

. :DRCUND(S*Constant*Range/ 1  2 

4.3) 

;  "KU" 

/  f*l  O  f  i 

*_  c,  0 

PRINT 

"AVER 

AGE  ENERGY* . 

.  .  " : DROUND*  Avarea*Constant*Ran 

as/ i 

24000 . 4 ) 

*’n  J" 

_  030 

PRINT 

"STANDARD  DEVIATION  OF 

TOTAL  ENERGY* . 

. . " : DROUND < Sar ea*Cons tan  t “Rang 

e/12 

4000  .?.);' 

•  >•« 

R  J 

_  0  4  0 

PRINT 

2050 

PRINT 

••**** 

•PHOTON  CONSTANT  = 

" : Constant ; " 

PYRO  CONSTANT'" ;Pyro_const ; "** 

..060 

PRINT 
DISP  1 
,T  r 

"TOTAL  Or  ":Bad:”  BAD 

PULSES  OUT  OF  " 

.•Pulse:  "RECORDED" 

.  080 

n  t  e  f  = 

701  THEN  Keys 

2  30 

PRINT 

CHRS  ( 

1 33) . "PROGRAM  NOW 

IN  GRAPH  MODE. 

USE  KEYS  TO  CONTINUE" .CHRSt 12 

[  N  T 


:  !  •  1  r  •!*!•!!!*  f  !!*!!!!!!!!! !  END  OF  -ULSE  C ALCULA7 IONS  !  1  !!!'!!  !!’!!’!!  ’  1  !  !  1  » 

I  I  »  »  t  I  »  t  f  I  »  *  »  t  »  I  I  I  I  t  I  t  t  I  I  >  t  «  t  I  >  II  I  t  I  II  t  t  t  I  t  l  t  I  *  »  t  1  t  I  I  I  t  I  I  I  I  I  !  I  I  I  I  t  I  I  !  '  «  •  I  t  >  t  t  t  t  I 


20 


20 

<  e  y  s  : 

j  *  j  0 

ON  <-  ■  :  _ 

^  !  4  0 

O'.  v  E 5  _ 

2150 

ON  KEY  E  _ 

w'  »  6  O 

u N  <L  " 

2 1  70 

ON  KE-  2  .. 

2180 

ON  KEV  *;  - 

2190 

ON  KE J  *  ■. 

2200 

ON  KE  -  r  . 

2210 

ON  KE  i  . 

2220 

ALPHA  I" 

2230 

!  LOOP  RE- 

2240 

GOTO  2220 

2250 

Screen:  ?*** 

2260 

GOSUB  : - 

2270 

graph:::  . 

2280 

INPUT 

2290 

GINIT 

2300 

GRAPHICS  : 

2310 

Plot:  ALPHA 

2320 

VIEWPORT  C 

2330 

Xmin=- .  5- 

2340 

Axis  i : p  =  I 

2350 

WINDOW  Xr,: 

2360 

CLIP  0.102 

2370 

AXES  100.2 

2380 

CLIP  OF" 

2390 

CSIZE  2.5. 

2400 

t  -.1 

2410 

LORG  6 

2420 

FOR  1=0  TO 

2430 

MOVE  1,0 

2440 

LABEL  I 

2450 

NEXT  I  ! 

2460 

•  Nil! 

2470 

LORG  8 

2480 

FOR  1=0  TO, 

2490 

MOVE  0.1 

2500 

LABEL  I 

2510 

NEXT  :  • 

2520 

MOVE  Inter 

2530 

LORG  6 

2540 

CSIZE  3.1 

2550 

LABEL  "TIM: 

2560 

MOVE  Xn in . i 

2570 

DEG 

2580 

LDIR  30 

2590 

LABEL  "FjWi 

_soo 

LDIR  0* 

26 1  0 

MOVE  0  .  C 

2  520 

FOR  1=’  'I 

2630 

draw  : 

*>  “  /.  .*• 

NEXT  * 

_  ^  5  0 

-CRG  5 

^  h,  J 

L  b  x  Z  c  2. .  i  . 

_  i. 

.-"OR  ,  =  ’  *" 

X  =  Rava-  :  . 

v.' 

Ffpw#  A-2.  Program  feting  (contimradl. 

.--'.EL  "VIEW  GRAPH"  GOTO  Screen 
'EfiD  PROGRAM"  GC'U  cimsh 
ARAMETERS"  GOTO  Params 
"CAMPLE  STATS"  GOTO  Statistics 
■  _JT  GRAPH"  GOTO  Hard 
"  OGGLE  PRil-iTER"  GOTO  Toggle 
"'"ORE  GRAPH"  GOTO  Disc 
"'JEW  SAMPLE"  GOTO  Start 
’.nANGE  PHOTON"  GOTO  Const 
"CHANGE  PYRO"  GOTO  Pyro 
:  'URNS  ON  ALPHA  NUMERIC  SCREEN 


ITS  PULSE  IN  KW  vs  nS  ON  CURRENT  PLOTTING  DEVICE  ********* 

!  SHUTS  OFP  USER  DEFINED  KEYS 
?  SHUTS  OFF  GRAPHICS  ! 

;E  NUMBER . Porapn 

!  INIt!AlIZE  GRAPHICS 

f  TURN  ON  GRAPHICS 

!  SHUT  OFF  PRINTOUT  (NOT  GRAPHICS) 

!  0.100  !  SET  LIMITS  OF  GRAPHING  AREA  (DISPLAY  UNITS) 

•ai*i 024  !  ALlOW  ROOM  FOR  LABELING  OF  Y  AXIS 
t ,al*Range*Constant/ 1280+1 >*10  !  ALTER  AXES  TO  SAMPLE 

cerval*1 024+25*Interval > , - . 1 *Ax is_top , Ax is_top  LSETS  RANGE 
?rval .O.Axis_top  !  DISABLES  GRAPHICS  OUTSIDE  RANGES 
.5.2  A  *  DRAWS  AXES 

!  ALLOWS  GRAPHICS  THROUGHOUT  VIEWPORT 
f  SETS  SIZE  OF  GRAPHICS  CHARACTERS 
ENG  OF  THE  X-AXIS _ 

!  SETS  ORIGIN  OF  GRAPHICS  CHARACTERS (TOP  CENTER) 

STEP  1000 

!  LABELS  TIME  AXIS  EVERY  1000  nS 


!  SETS  ORIGIN  OF  GRAPHICS  CHARACTERS  (RIGHT  CENTER) 
Axis_top  STEP  50 

!  LABELS  INTENSITY  AXIS  EVERY  50  KW 


vai*500.-'.05*Axis_toD  !  X  AXIS  LABEL  LOCATION 

!  SETS  ORIGN  OF  GRAPHICS  CHARACTERS  (TOP  CENTER) 

!  SETS  SIZE  OF  GRAPHICS  CHARACTERS 

“  ( "i  S ) 11 

Sxisltop/2  !  Y  AXIS  LABEL  LOCATION 

!  ROTATE  GRAPHICS  CHARACTERS  90  DEGREES 

E*  <  KW)" 

?  RETURN  GRAPHICS  CHARACTERS  TO  HORIZONTAL 
!  PLOTING  PULSE 


'SET  GRAPHICS  CHARACTER  ORIGIN  (CENTER) 
! SET  SIZE  OF  GRAPHICS  CHARACTERS 
•MARKING  SIGNIFICANT  POINTS  WITH  "*" 

•  ::aoh) 

:j  THEN  27 i If!  !  IF  OUT  OF  BOUNDS  THEN  DON’T  MARK 


Copy  crmflaWe  to  DTIC  doe*  not 
permit  fully  legible  reproduction 


2780 
2730 
2800 
2810 
2820 
2830 
2840 
2850 
2360 
2870 
2880 
2890 
2900 
2310 
2920 
2930 
2940 
2945 
f  2950 

2960 

2961 

2962 
-970 
2980 
2990 
3000 
2010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
2160 
3170 
3180 
3190 
8200 
3210 
3220 
3230 
3240 
3250 


►  ft*’-#***********-*-*'* 


Figure  A-2.  Program  lifting  (eontinuod). 

Y* A(  Pg rapn  .  X  > *Ran*pe*Cons  t an  t /I  28  !  CONVEr  TO  kIlOwH*  "Ls 

IF  Y<G  OR  Y>Axis_top  THEN  2740  *  REJECT  IF  OFF  GRAPH 

MOVE  X*Intervai , Y 
LABEL  "*" 

NEXT  J 
MOVE  0.0 

PLOTTER  IS  3. "INTERNAL"  !  SETS  PLOTTER  TO  CRT  IN  CASE  EXTERNAL 

GOTO  Keys  1  ■•****»»****•** *»■»»*««**  END  OF  GRAPH  SECTION  . . 

**#**'****'*-***-*-**r4*  +  *i**’***-*-*i****‘#***  +  ***i**#**lf4Hf+’+  +  +  *4Hf4Hf-*it-innn>**irlnnt*i*-l**-0’-+*it-»+- 

Hard:  PLOTTER  IS  705."HPGL"  ! »*»*****»*. ENABLES  EXTERNAL  PLOTTER—*  —  * 

GOSUB  Off  keys 

DISP  "PREPARE  PLOTTER  AND  CONTINUE” 

PAUSE 

GOTO  Plot  ?****»»»*«»»******»*i»***»***»«**««**^< **•****«*********•**«*»*< 
Toggle:  !***.*»**  TOGGLES  PRINT  DEVICE  ( INTERNAL/EXTERNAL)— *************** 
GOSUB  Of f_keys 
IF  Pr inter=701  THEN  2900 
PRINTER  IS  701 
Pr  inter=701 
DISP  "PRINTER  IS  70tM 
GOTO  2930 
Pr inter= 1 
PRINTER  IS  1 
DISP  "PRINTER  IS  1" 

GOTO  Keys 

Params:’— — — *  PRINTS  OUT  SAMPLE  STATS  WITH  CURRENT  CONSTANTS***** 
GOSUB  Of f_k eys 

INPUT  "INPUT  THE  PULSE  YOU  WANT  TO  SEE  STATS  0N".P 
IP  P<=Pulse  AND  P>0  THEN  2970 

DISP  "SORRY. . .ILLEGAL  PULSE  NUMBER.  MUST  BE  BETWEEN  0  AND  ";Puise 
GOTO  3050 

Tstart=Parameter< 1 .P)  f  SETS  PARAMETERS  FROM 

Tmax=Parameter <2 ■ P >  !  PARAMETER  ARRAY 

Tedge=Parameter(3 ,P> 

Tend ’Parameter <4 ,P) 

T 1 ’Parameter <5 ,P> 

T2=Parameter <6 .P) 

Max=A<  P . Tmax ) 

GOSUB  Printout 

GOTO  Keys  {»»******»**«»—»«»«*****«*«•*•<»•*•**»,*«**«**«■*•*<****««■*«***•*•»■* 
Disc:  ;**,**«*•**»*****. STORES  GRAPH  ON  DISC  ********•..-•**.**.*•.***•.**«. 

GOSUB  Of  f_keys 

INPUT  "ENTER  THE  FILE  TO  STORE  PULSE  0N".DS 

INPUT  "WHAT  PULSE  DO  YOU  WANT  TO  STORE?" .$ tor e_puise 

MASS  STORAGE  IS  ": INTERNAL , 4 , 1"  »  MASS  STORAGE  DEVICE  LEFT  DISC  DRIVE 

CREATE  BOAT  DS.2051,8  !  CREATE  DATA  FIlE  O'7  APROPRIATE  LENGTH 

ASSIGN  'sPath  TO  DS  !  ASSIGN  OUTPUT  PATH  TO  DATA  FILE 

FOR  1=1  TO  N 

Tempy ( I ) =A( Store_pulse . I > *Cons tant’Range/ 1 28  IDEFIME  TEMPORARY  ARRAY  <KW> 
NEXT  I 

OUTPUT  ^Path ; 1 020  !  STORE  NUMBER  UE  POINTS  FOR  GRAPHICS  PROGRAM 

OUTPUT  «Path: Interval  !  STORE  nS  PER  ARRAY  ELEMENT 
OUTPUT  iPath : Tempy< *>  •  STORE  ARRAY  <KW>  UNFORMATTED 
MASS  STORAGE  IS  " : INTERNAL"  »  MASS  STORAGE  DEVICE  RIGHT  DISC  DRIVE 
PRINT  "PULSE  STORED  <KW  vs.  nS>  ON  FILE  NAMED  ":DS 

GOTO  Keys  . . . . . . 

Const :  f ***«*»»***«*»**»»»»CHANGES  PHOTON  DRAG  CONSTANT*  — —  —  —  *  — *  — 

GOSUB  Of f _k eys 

DISP  "USE  KEYS  TO  CHOOSE  METHOD  OF  ENERGY  DE TERM 1NAT IrN" 

ON  KEY  1  LABEL  "CALORIMFTER"  GCTC  328.;, 

22  ✓ 


F%uf«A-2.  Profram  lifting  (oont>nu«d). 


ON  KEY  2  LABEL  "PYRO  ELECT."  GOTO  3330 

GOTO  2270 

.2  S'.-  GO  SUB  Of t  _<  eys 

2230  INPUT  "ENTER  AVERAGE  CALORIMETER  READING  IN  mJ“.  M i i i joules 

Constant  =Const an t*Mi  i  i  joul  es/<  Avarea*Cons t ant "Range /I 2 AO 00  > 

;2U>  PRINT  "NEW  PHOTON  DRAG  CONSTANT-" :DROUND<Constant .5) 

2 3 2 1 1  GOTO  Keys 

:-33U  GO  SUB  Off  Leys 

33-0  INPUT  "HOW  MANY  ENERGIES  DO  YOU  WANT  TO  TAKE"  .Pyro_rei- 

3350  DISP  "CHANGE  INPUT  TO  PYRO  AND  CHANGE  TIMEBASE  TO  FIND  MAXIMUM. THEN  <C0NT 

;  NuE 


3360 
2370 
3380 
3390 
3  A00 
3  A !  0 
3A20 

4  3  '! 

3  A  AO 
7A50 
3460 
;  a  7  0 
2  a  SO 
•  a9U 
3500 
25 ;  0 
3520 
2530 
3540 


PAUSE 

DISP  "READY  FOR  PULSE  FROM  PYRO" 

Ave_py ro=G  !  INITIALIZE  AVERAGE  PYRO  READING 
FOR  J=1  TO  Pyro_rep 

OUTPUT  704:8  •  ARMS  THE  6500  DIGITIZER 

OUTPUT  704; 11  !  PUTS  THE  6500  INTO  OUTPUT  MODE 

ENTER  704  USING  !  ENTER  BUFFER  ARRAY  UNFORMATTED 

OUTPUT  704 ; 0  !  CLEARS  THE  6500  OF  DIGITAL  OUTPUT 

DISP  J:"  OF  "  ;Py ro__rep  : "  SHOTS" 

Max  _pyro=0 
FOR  1=1  TO  500 


IP  B < I ) <  =  !  2 7  THEN  3490 
Bi I )-B<  D-256 
B( I ) =B< I ) -Dc_of  f  set 
IF  Max_pyro>B< I >  THEN  3520 
Max_py  ro=B< I ) 

NEXT  I 


*  CONVERT  TO  WORD  2'S  COMPLEMENT 

*  FROM  BYTE  2'S  COMPLEMENT 
!  ACCOUNT  FOR  DC  OFFSET 

*  FINDS  THE  MAXIMUM  VOLTAGE 


Ave_pyro=Max_pyro/Pyro_rep*Ave_pyro  f  FINDS  THE  AVERAGE  MAXIMUM  VOLTAGE 

NtXT  j 


2560 

2570 

3580 


Znergy=Ave_pyro*Rar.ge*Pyro_const/128  !  FINDS  THE  ENERGY  FROM  PYRO  CONST 
Constant»DROIJND(tnergy/ (Avar ea»Range/ 128000) ,6)  ?NEW  PHOTON  DRAG  CONSTANT 
PRINT  "NEW  CONSTANT-"; Constant 

GOTO  Keys  !  **»*»»*»*»■»»«•«»«*«»»*««**»*»»»****»»««**•>*««■•***»■»*»**»»»»*»»* 
Offjceys:  !  SHUTS  OFF  USER  DEFINED  KEYS*‘»-**»**»****»*»**»**»»*»»»»«««*«** 

OTHERWISE  PROGRAM  MAY  BRANCH  IN  MIDDLE  OF  SUBROUTINE 


3600 

OFF 

KEY 

0 

26 1  0 

OFF 

KEY 

1 

3620 

OFF 

KEY 

2 

2630 

OFF 

KEY 

3 

3640 

OFP 

KEY 

4 

:  650 

OFF 

KEY 

5 

3660 

OFF 

KEY 

6 

2670 

OFF 

KEY 

7 

ipyo 

OFF 

KEY 

8 

•  690 

OFF 

KEY 

9 

3  7  0  0 

RETURN 

**•*■*•*« 

3710  Pyro:  !  DETERMINES  NEW  PYROELECTRIC  CUNSTANT***»»*»» »****« •*****»*»»**• 

1720  GC'SUB  Off  keys 

3730  PRINT  "OLD  PYRO  CONSTANT-" :Pyro_const 

:740  PRINT  "ENTER  PYROELECTRIC  PULSE  THRU  DIGITIZER" 

3750  INPUT  "HOW  MANY  SHOTS  DO  YOU  WANT  TO  TAKE”. Shots 

3760  Avmaxvol t*0 

3770  DISP  "READY  FOR  DIGITIZER  INPUT.  HIT  (CONTINUE)  WHEN  READY" 

3  T  PC  PAUSE 

:7S:  FOR  J=1  TO  Shots 

7:  OUTPUT  704:8  !  ARMS  THE  6500  DIGITIZER 

■  •  OUTPUT  704:11  !  PUTS  THE  6500  INTO  OUTPUT  MODE 

ENTER  704  USING  !  ENTERS  BUFFER  ARRAY  FROM  6500 

OUTPUT  704 ;0  *  CLEARS  THE  6500  OF  OUTPUT  MODE 
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Figure  A-2.  Program  Meting  (continued). 


8840 

:350 

i860 

3870 

'.880 

3890 

3900 

3910 

8920 

3930 

3940 

39S0 

es 

8960 

3970 

3980 

?990 

4000 

^010 

4020 

4030 

4040 

j041 

4050 

“060 

4070 

4080 

4090 


Haxvo 1 1  =  0 
F'JR  1=1  T 1 J  j00 
IF  B< I ) <= 1 27  THEN  8880 
B< I)-B<I>-256 
B(I>*B<D -Dc_of  f  set 
IF  BdXMaxvoit  THEN  39"  0 
Max vo 1 t  =  B(  I ) 

NEXT  I 

DISP  J : "OF" : Shots ; "SHUTS  T2  EE  TAKEN" 

Avmaxvo i  t “Maxvo  1  t/Sho ts+A vr.ax vo i  t  !  FINDS  AVERAGE  MAX  PYRO  VOLTAGE 


! CONVERT  73  WORD  2'S  COMPLEMENT 
f  "ROM  BYTE  2'S  COMPLEMENT 

!  ACCOUNT  FOR  DC  OFFSET 
•  FINDS  MAXIMUM  PYR0  VOlTAGE 


NEXT  J 

INPUT  "ENTER  ENERGY  READING  FROM  CALORIMETER  IN  MILIJOULES  ".Milijoul 


Pyro  const*DROUND<Mi 1 1 iou!es/< Avmaxvol t*Ranoe/ ! 2S) .5) 

PRINT  "NEW  PYRO  CONSTANT-" ;^yro_const  •  HEW  PYRO  CONST.  CALIBRATED  TO 

CALORIMETER  ■ 

ON  ERROR  GOTO  4010  1 

Constant”DROUND(Mi 1 i joules/ < Avarea*Range/ 1 28000 ). 6 >1  NEW  PHOTGN  DRAG  CONST 
GOTO  4020 

PRINT  "NO  ENERGY  READINGS  AVAILABLE  TO  CHANGE  PHOTON  DRAG  CONSTANT" 

OFF  ERROR 

GOTO  Keys! ******************** **«.*»***«* . . 

Finish:  PRINTER  IS  1 

PRINT  CHRS< 12) 

GRAPHICS  OFF 

DISP  "HIT  <CONTINUF>  TO  LOAD  PROGRAM  DIRECTORY" 

PAUSE 

LOAD  "AUTOST" 

END 
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Y  OUTPUTS 
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1026 

Data  bits 
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Out 
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Input  RanRe 
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*7 
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Low: DC  High:AC  Low:+  High:- 

Low : INT  High: EXT 

*6 

0.25 

64 

MS 

*5 

♦ 

0.50 

32 

MS 

*4 

} 

1.0 

16 

NS 

y3 

2.0 

8 

8 

*2 

LSB 

5.0 

4 

4 

*1 

LOW 

LOW 

2 

2 

*0 

LOW 

LOW 

1 

1 

Logic 

POS.  TRUE 

NEG.  TRUE 

POS.  TRUE 

NEG.  TRUE 

Code 

2's  Comp 

1  of  5 

1>  Pull  Scale 

Weighing: 

Defin. 

AC/DC 

Signed  Binary 

0=2  4=50 

1=5  5=100 
2=10  6=200 
3=20  7=500 

8=1000 


Figure  A-4.  Raadingi  from  digitizer. 


c.  Input  Range  (lines  360  to  400).  In  this  case,  only  bits  2  to  6  are  of  interest. 
Since  these  hits  are  negative  true  logic,  the  low  bit  corresponds  to  the  input  range.  Line  320 
calculates  the  input  range  for  Y3  to  Y6.  The  exception  to  the  equation  (Y2)  is  trapped  in  line 
400. 


d.  Offset.  The  offset  is  calculated  in  line  360.  The  binary  AND  of  the  offset  byte  with 
127  masks  out  the  sign  bit  to  give  the  absolute  offset  in  percent.  Divide  by  100  to  get  the  frac¬ 
tion  of  input  range.  Multiply  by  —  1  A  (Y7)  to  account  for  sign. 

e.  Interval  (lines  290  to  350).  Bits  Y„  to  Y3  are  of  interest  here.  For  this  applica¬ 
tion,  internal  time  base  in  nanoseconds  is  assumed.  Line  270  determines  the  weighing  number. 
The  binary  AND  of  15  and  the  interval  byte  gives  only  Y„  to  Y3.  The  EX-OR  of  15  and  Y0  to 
Y3  gives  the  desired  positive  true  weighing  factor.  Since  there  is  no  simple  equation  to  generate 
the  interval  weight,  a  loop  is  used  to  read  the  proper  weighing  factor.  An  array  could  be  used, 
but  would  require  more  memory. 

f.  Offset  Check  (lines  540  to  610).  There  is  no  simple  way  to  calculate  the  step 
offset  from  the  percent  offset  to  better  than  an  interval  of  4  (percent).  To  eliminate  this  possi¬ 
ble  3-percent  error  in  voltage,  the  program  prompts  the  user  to  arm  and  trigger  the  digitizer 
with  O  V  input  to  check  the  calculated  step  offset  with  the  true  screen  offset.  If  they  differ  by 
more  than  4  (percent),  the  program  re-prompts  the  user.  If  they  are  within  4  (percent),  the  true 
screen  offset  is  used  in  the  program. 

3.  Array  Definitions  and  Pulse  Data  Entry  (lines  700  to  890). 

a.  Integer  Arrays  A(120,  1023).  Holds  pulses  in  interval  step  units. 

•  Room  for  120  pulses  with  1024  samples  each. 

•  B(1026).  Buffer  array  for  fast  data  entry  by  unformatted  ENTER  statement. 

•  Parameter  (6,120).  Holds  the  times  of  the  six  important  events  of  each  of  120 
pulses. 

•  Bad  pulse  (120).  Marks  the  bad  pulses  to  be  excluded  from  statistical 
analysis. 

b.  Real  Arrays  Area  (120).  Area  under  pulse  in  interval  step  units. 

•  Aspike  (120).  Area  under  spike  in  interval  step  units. 
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•  Tempy  (1023).  Temporary  storage  for  pulse  intensity  in  kW  used  for  pulse 
disc  storage. 

NOTE:  Most  calculations  and  storage  are  left  in  interval  step  units  to  reduce  run 
time  and  memory  use.  Conversion  to  proper  units  occurs  only  for  I/O. 

•c.  Pulse  Data  Entry.  The  HP  9836  computer  uses  an  HPIB  to  interface  with  the 
6500  digitizer.  The  digitizer  address  in  this  program  is  704.  There  are  four  digitizer  I/O  com¬ 
mands  necessary  for  this  program. 

•  OUTPUT  704;8.  Arms  the  6500  Digitizer  (single  sweep  mode). 

•  OUTPUT  704;11.  Puts  the  6500  in  output  mode. 

The  6500  now  waits  for  a  trigger  pulse  and  then  records  waveform  data. 

•  ENTER  704  USING  Enters  the  data  from  the  6500  into  a  buffer 

array  (unformatted).  This  buffer  is  necessary  because  an  unformatted  read  into  a  two- 
dimensional  array  is  not  allowed,  and  single  element  entry  into  a  two-dimensional  array  takes 
much  longer. 

•  OUTPUT  704;0.  Clears  the  6500  of  digital  output. 

The  buffer  array  is  now  transferred  into  the  two-dimensional  A (120,  1023)  array. 

NOTE:  Data  stored  in  the  A(120,1023)  array  is  in  2's  complement  bytes.  No  fur¬ 
ther  processing  is  done  at  this  time,  because  pulse  rate  is  already  limited  to  0.7  Hz. 

4.  Convert  to  Word  2’s  Complement  and  Subtract  DC  Offset  (lines  1000  to 
1080).  This  program  segment  converts  the  byte  2’s  complement  data  from  the  digitizer  to  the 
word  2's  complement  used  by  the  computer.  The  DC  offset  is  subtracted  from  each  pixel  (see 
Block  2).  The  pulse  is  also  rectified  by  taking  the  absolute  value,  so  that  either  inverting  or 
non-inverting  amplifiers  can  be  used  after  the  detector. 

5.  Find  Maximum  of  all  Pulses  (lines  1050  to  1060).  The  maximum  intensity  of  the 
sample  of  pulses  is  found,  so  that  any  pulses  below  40  percent  of  this  value  can  be  rejected. 
This  maximum  is  also  used  to  determine  the  intensity  range  for  graphing  the  sample  pulses. 

6.  Find  Pulse  Start  Time  and  Maximum  Value  (lines  1160  to  1240).  Ideally,  the 
start  of  the  pulse  should  be  recognized  by^a  change  in  slope.  However,  noise  frequency  limits 
this  option.  The  method  used  in  the  program  triggers  from  zero.  The  program  keeps  track  of 
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where  the  waveform  last  hit  zero  (line  1200).  If  the  waveform  does  not  return  to  zero  after  200 
sample  units,  the  last  zero  is  considered  the  pulse  start  (see  Figure  A-5).  This  step  of  200  units 
may  have  to  be  decreased  for  short  pulses  or  long  time  bases.  The  program,  then,  records  the 
maximum  intensity  (Max)  and  the  corresponding  time  Tmji. 

7. .Forty  Percent  Reject  Level  (line*  1250  to  1290.)  If  a  pulse  maximum  (Max)  is 
less  than  40  percent  of  the  sample  maximum  (Maxval),  then  the  pulse  is  called  “bad”  and  is 
not  included  in  statistical  pulse  analysis.  This  rejection  eliminates  false  triggers  caused  by  ex¬ 
ternal  corona  discharge  or  turning  detectors  on  or  off. 

8.  Locate  Spike  End  From  Slope  Change  (lines  1320  to  1370).  The  spike  and 
(Tedge)  is  defined  at  the  point  where  the  slope  of  the  waveform  suddenly  becomes  less  negative 
and  rises  above  a  certain  value.  This  equation  is  in  line  1360  and  may  have  to  be  tailored  for 
different  time  bases.  A  pulse  can  be  rejected  at  this  point,  if  the  timing  sequence  is  incorrect 
(line  1400). 

9.  Locate  Pulse  End  (lines  1400  to  1440).  The  pulse  end  is  defined  as  the  first  point 
after  the  spike  end  where  the  waveform  drops  to  zero  T^  .  If  the  waveform  does  not  drop  to 
zero,  the  pulse  is  rejected. 

10.  Find  Spike  Energies  and  Hsdf  Widths  (lines  1470  to  1560). 


Since  the  intensity  is  a  stepped  function,  rectangular  integration  is  used  to  determine 
pulse  energy.  This  segment,  also,  finds  T,  and  T2  to  determine  the  pulse  half  width  (Figure 
A-5). 

11.  Significant  Parameters  Stored  (lines  1720  to  1790).  The  times  at  which  the  six 
significant  events  occur  for  each  pulse  are  stored  in  the  Parameter  (6,120)  array.  These  times 
are  used  later  to  find  corresponding  intensities  and  temporal  characteristics. 

12.  Print  Significant  Parameters  (lines  1590  to  1710).  Prints  single  pulse  statistics 
to  current  printing  device.  Converts  to  kW,  mj,  and  ns  on  printout  using  current  photon  drag 
constant. 

13.  Calculate  Sample  Statistics  (lines  1810  to  1970).  This  program  section  looks  at 
all  pulses  entered.  If  the  pulse  is  bad  because  of  low  intensity  or  a  timing  error,  it  is  not  includ¬ 
ed  in  the  analysis.  Energy  and  maximum  intensity  averages,  and  standard  deviations  are 
calculated  for  the  group  of  pulses. 


14.  Key  X  (lines  2120  to  2240).  At  this  point  all  major  calculations  are  completed.  The 
photon  drag  or  pyroelectric  constants  can  be  changed,  or  data  can  be  sent  to  external  output 
devices.  By  hitting  a  user-defined  key,  the  program  executes  a  subprogram,  returns  to  this  spot, 
and  waits  for  another  key  to  be  depressed.  Keys  are  disabled  during  subprograms  so  output 
can  be  completed. 

15.  Key  0  “VIEW  GRAPH”  (lines  2250  to  2770).  This  key  produces  a  graph  of  power 
vs.  time  of  the  desired  pulse  number  on  the  computer  CRT.  The  top  of  the  power  axis  is  deter¬ 
mined  by  the  maximum  sample  intensity.  The  significant  points  of  the  laser  pulse  are  marked 
with  a  for  reference.  The  program  then  returns  to  Key  X. 

16.  Key  1  “PLOT  GRAPH”  (lines  2780  to  2820,  2310  to  2770).  This  key  transfers 
control  to  an  external  plotter  and,  then,  branches  into  Key  0  to  produce  a  plot  of  power  vs.  time 
on  paper. 

1 7.  Key  2  “TOGGLE  PRINTER”  (lines  2830  to  2930).  This  key  is  a  fast,  one  step 
way  of  changing  the  printer  from  internal  to  external  or  vice  versa  (i.e.,  CRT  display  or  hard 
copy).  Otherwise,  the  program  would  have  to  be  Btopped,  a  command  executed  and,  then,  the 
program  continued. 

18.  Key  3  “STORE  GRAPH”  (lines  3060  to  3210).  This  key  stores  on  disc  the  pulse 
intensity  in  kW.  The  number  of  points  and  the  number  of  nanoseconds  per  array  index  is  also 
stored.  The  format  of  the  data  file  is  compatible  with  a  graphics  program  to  make  presentation 
quality  graphs. 

19.  Key  4  “NEW  SAMPLE”  (line  2190.)  This  key  prepares  for  a  new  pulse  sample  but 
does  not  reinitialize  the  photon  drag  and  pyroelectric  constants.  The  digitizer  control  settings 
can  be  changed  for  the  new  sample.  Program  continues  on  line  80. 

20.  Key  5  “PARAMETERS”  (lines  2940  to  3050).  This  key  initializes  the  time 
markers  to  those  of  the  desired  pulse  from  the  parameter  (6,120)  array.  The  program  then 
prints  out  the  pulse  rise  time,  half  width,  spike  duration,  total  duration,  spike  energy,  total 
energy,  and  peak  intensity  using  the  current  photon  drag  constant.  This  printout  is  sent  to  the 
current  printing  device  as  specified  by  Key  2. 

21.  Key  6  “SAMPLE  STATS”  (lines  1980  to  2070).  This  key  prints  out  the  pulse  sam¬ 
ple  averages,  and  standard  deviations  of  pulse  energy  and  maximum  intensity  on  the  current 
printing  device.  The  energies  and  intensities  are  computed  using  the  current  photon  drag  and 
pyroelectric  constants  which  are,  also,  printed  out. 
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22.  Key  7  “CHANGE  PHOTON'”  (lines  3220  to  3580|.  This  key  is  used  to  change  the 
photon  drag  constant.  If  the  average  energy  computed  by  the  program  does  not  match  a 
calorimeter  reading,  the  photon  drag  detector  constant  can  be  changed  by  entering  a 
calorimeter  reading  or  by  using  the  pyroelectric  detector  as  a  calorimeter.  The  photon  drag  con¬ 
stant  is  in  units  of  kW/V. 

23.  Key  8  “CHANGE  PYRO”  (lines  3710  to  4030}.  This  key  calculates  a  new 
pyroelectric  constant  in  tnJ/V.  The  average  maximum  voltage  from  the  pyroelectric  detector 
through  the  6500  digitizer  is  determined.  This  average  voltage  is  proportional  to  the  average 
energy  which  is  entered  in  mj  to  determine  the  constant  in  mJ/V  of  the  pyroelectric  detector. 
The  pyroelectric  detector  can  now  be  used  as  a  standard  to  calibrate  the  photon  drag  detector. 

24.  Key  9  “END  PROGRAM”  (lines  4040  to  4090).  This  key  reloads  the  disc  pro¬ 
gram  directory. 

General  Comments.  Some  of  the  algorithms  used  in  this  program  sacrifice  speed  in 
order  to  reduce  memory  requirements.  In  the  current  program,  pulse  analysis  is  limited  to  1 20 
pulses.  This  limit  could  be  considerably  increased  if  comments  or  unneeded  segments  are 
eliminated.  LOAD  SUB  and  PURGE  SUB  commands  could,  also,  be  used  to  load  and  purge 
program  segments  as  they  are  needed  on  the  run. 

For  lifetime  tests,  a  loop  could  be  put  in  Section  3  (Pulse  Data  Entry)  to  record  every 
10th  or  100th  pulse  to  increase  sample  range. 

All  pulse  arrays  in  this  program  are  kept  in  step  units  until  the  time  of  output.  This  is 
done  for  several  reasons.  The  array  as  integer  step  units  takes  up  less  room  than  the  real  array 
converted  to  kW.  The  processing  of  integer  numbers,  also,  takes  less  time  than  floating  point 
operations.  Another  advantage  is  that  when  the  photon  drag  constant  is  changed,  none  of  the 
pulse  storage  arrays  have  to  be  accessed  since  conversion  is  done  at  output. 

In  some  cases,  program  efficiency  is  given  up  for  clarity.  After  the  program  is 
understood,  efficiency  can  be  improved  by  removing  comments  or  I/O  prompts  and  combining 
lines. 
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